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^7  This  paper  presents  a  nontechnical  review  of  current  research  efforts  in 
the  mathematical  theory  of  viscoelastic  materials.  Recent  results  concerning 
the  existence,  uniqueness  and  regularity  of  solutions  for  initial  value 
problem  as  well  as  for  steady  flow  problem  are  discussed  and  a  number  of 
open  problems  are  pointed  out.  0  ^  ^  Atfyi-cWs  ■  r'Mvd  fi  ; 


AMS  (MOS)  Subject  classifications!  35K 15,  35K20,  35K55,  35L15,  35L67, 

35L70,  35M05,  35S10,  45K05,  76A10 


.L 


Key  Words:  ^Viscoelastic  fluids,  initial  value  problems-^  steady  flows, 
nonlinear  evolution  problems,  energy  estimates,  development 
of  shocks , ^singular  kernels  .  _ 

Work  unit  Numbers  1  (Applied  Analysis) ,  2  (Physical  Mathematics) 


SIGNIFICANCE  AND  EXPLANATION 


This  paper  is  based  on  a  lecture  presented  at  the  Symposium  on 
Viscoelasticity  and  Bheology,  Madison/  October  1984.  The  lecture  was  aimed  at 
explaining  to  a  general  audience  some  of  the  mathematical  problems  associated 
with  model  equations  for  viscoelastic  materials  and  giving  an  impression  of 
the  current  state  of  efforts  in  this  field.  Questions  of  existence, 
uniqueness  and  regularity  of  solutions  to  initial  value  problems  and  to  steady 
flow  problems  are  discussed  and  some  open  problems,  which  the  author  considers 
to  be  of  particular  interest,  are  pointed  out. 
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RECENT  DEVELOPMENTS  AND  OPEN  PROBLEMS  IN  THE 
MATHEMATICAL  THEORY  OF  VISCOELASTICITY 


Michael  Renardy 

1.  EXISTENCE  RESULTS  FOR  INITIAL  VALUE  PROBLEMS 

The  existence  of  solutions  for  systems  of  ordinary 
differential  equations 

ft  «  f(u),  u< 0)  -  u0  ,  (1.1) 

can  be  established  by  using  the  iteration 

0n+1  -  f(un),  un+1(0)  *  u0  ,  (1.2) 

which  can  be  shown  to  converge  for  small  enough  t  >  0.  For 
partial  differential  equations  such  as 

ut  ■  ♦<UX)X'  u(t  “  °'x>  “  uq( x)  ,  (1.3) 

the  iteration  (1.2)  cannot  be  used.  This  is  because  at  each 
step  of  the  iteration  we  would  be  talcing  two  x-derivatives 
and,  if  data  are  sufficiently  oscillatory,  these 
x-derivatives  will  grow  larger  and  larger.  An  iteration 
that  can  be  used  is  the  following 

uj+1  *  ♦’  u”+1<t  *  0,x)  »  uQ(x)  .  (1.4) 

This  example  illustrates  a  general  strategy  used  for  so- 
called  "quasilinear"  equations.  In  such  equations  the 
highest  derivatives  occur  only  linearly,  and  nonlinearities 
contain  only  lower-order  derivatives.  One  then  puts  an  n  + 
1  on  the  highest  derivatives  (the  "principal  terms”  in  the 
equation)  and  an  n  on  lower  order  derivatives.  Two  things 
must  then  be  proved:  first  that  the  linear  problems  to  be 
solved  at  each  step  are  well-posed  and  secondly,  that  the 
iteration  converges. 


In  order  to  apply  these  ideas  to  problems  in  visco¬ 
elasticity,  we  have  to  know  that  the  equations  can  be 
written  in  some  quasilinear  form  and  we  have  to  know  what 
the  principal  terms  are.  In  the  following,  we  illustrate 
this  for  a  class  of  one-dimensional  model  equations,  whose 
structure  is  typical  of  problems  arising  in  viscoelasticity. 
Let  us  consider  equations  of  the  form 


utt  ■  “£<ux'uxt>x  +  9(ux>. 


+  /  m(t  -  T)h(u  (t),u  (t))  dT 


(1.5) 


with  appropriately  smooth  boundary  and  initial  data.  E'or 
m  >  0,  the  term  of  highest  differential  order  on  the  right 
hand  side  is  the  one  involving  uxxt.  If  the  derivative 
of  f  with  respect  to  the  second  argument  is  positive, 

(1.5)  can  be  regarded  as  a  parabolic  equation  and  the 
iteration  to  be  used  is 


“St1-  “D2£(uS-uSt,u«t  ♦  “Dif(u!'uSt)uSx 


+  g(ui)x  +  /  *<t  -  T>h(uJ(t),u"(T))xdT  .  (1.6) 


Under  appropriate  smoothness  assumptions  it  can  be  shown 
that  such  a  scheme  is  well-defined  and  converges.  This  is 
accomplished  by  adapting  a  result  due  to  Sobolevskii  [36] 
for  equations  without  integral  terms.  At  each  step  of  the 
iteration,  we  have  to  solve  a  linear  parabolic  equation  with 
time-dependent  coefficients.  In  Sobolevskii 's  approach, 
this  equation  is  regarded  as  a  perturbation  of  the  problem 
with  time-independent  coefficients,  evaluated  at  the  initial 
time  t  *  0,  and  a  convergent  iteration  based  on  the 
variation  of  constants  formula  can  be  formulated.  Finally, 
it  is  proved  that  the  mapping  un  ♦  un+*  given  by  (1.6)  is 
a  contraction  in  an  appropriate  Banach  space. 


One-dimensional  model  problems  of  parabolic  type  are 
discussed  e.g.  in  [11]/  [29].  In  [30]/  the  same  ideas  are 
extended  to  three-dimensional  flows  of  incompressible  fluids 
under  both  displacement  and  traction  boundary  conditions. 

If  u  *  0/  m,  g  and  h  are  smooth,  and  m,  g*  and 
Djh  are  positive,  then  (1.5)  can  be  regarded  as  a 
hyperbolic  equation.  The  appropriate  iteration  method  is 


n+1 


‘tt 


9*  (ux)uxx1+  /  m(t  “  OD^tuJJm/uJtTnuJJ^tt) 

—OP 


+  D2h(u5J(t),uJ(T))uJx(T)dT  .  (1.7) 

It  looks  surprising  at  first  that  we  can  put  an  n  on  the 
uxx(T)-term,  since  this  term  appears  to  be  of  the  same 
differential  order  as  the  uxx(t)-term.  However,  when  we 
differentiate  the  equation  with  respect  to  t,  then  uxx(t) 
becomes  uxxt,  but  uxx(t)  remains  uxx(t).  This 
explains,  on  a  heuristic  level,  why  this  term  can  be 
regarded  as  being  of  lower  order. 

At  each  step  of  the  iteration,  one  now  has  to  solve  a 
linear  hyperbolic  equation  with  time-dependent  coefficients. 
Existence  of  solutions  for  such  problems  can  be  shown  e.g. 
by  the  implicit  Euler  scheme.  Again,  a  contraction  argument 
is  used  to  show  convergence  of  the  iteration.  This  argument 
consists  essentially  of  two  steps:  First,  it  is  shown  that 
the  iterates  remain  bounded  in  a  certain  norm  (e.g.  that  all 
third  derivatives  of  un  have  bounds  independent  of  n). 
Then  one  uses  this  fact  to  show  that  the  iterates  converge 
in  a  weaker  norm  (i.e.  one  proves  that  second,  rather  than 
third,  derivaties  of  un  converge).  An  abstract 
formulation  of  these  ideas  is  given  in  the  work  of  Kato 
[17],  [19],  [20]. 

Applications  to  problems  in  viscoelasticity  are  due  to 
Mac  Camy  [23],  Staffans  [37],  Dafermos  and  Nohel  [5],  [6], 
Hrusa  [13],  Kim  [21],  Heard  [12]  and  the  author  [31],  [32]. 
Three-dimensional  problems  with  incompressibility  have  been 


studied  121] ,  [32] ,  but  only  for  the  pure  Cauchy  problem 
poseo  on  all  of  space,  not  for  initial-boundary  value 
problems . 

There  is  some  interest  in  relaxing  the  assumption  that 
m  is  smooth.  The  first  term  on  the  right  hand  side  of 
(1.5)  can  be  thought  of  as  arising  from  setting  m  equal  to 
the  derivative  of  the  6-function.  It  is  only  reasonable  to 
think  that  m  may  also  have  weaker  singularities  as 
t  —  t  ♦  0.  If  m  is  unbounded  at  0,  the  iteration  (1.7) 
will  not  work.  Instead,  one  will  have  to  put  n  +  1  on  the 
uxx( T)-term  as  well.  This  means  that  the  "quasilinearized" 
problems  to  be  solved  at  each  iteration  are  now  history 
dependent.  (In  both  the  parabolic  and  hyperbolic  cases  the 
linear  problems  involved  no  history  dependence  at  all,  i.e. 
history  dependence  was  basically  treated  as  a  perturbation.) 
In  a  recent  paper  [15J,  W.  Hrusa  and  the  author  show  that 
such  an  iteration  converges  tor  the  special  problem 

t 

utt  *  ♦(ux)x  *  1  m(t  ”  t)*(ux(t) )xdT  ,  (1.8) 

m 

where  +',  and  j  m(t)dT  are  positive,  m  is 

0 

positive  and  monotone  decreasing  and  has  an  integrable 
singularity  at  0.  The  method  used  is  based  on 
approximating  the  singular  kernel  by  smooth  kernels  and 
showing  that  it  is  possible  to  pass  to  the  limit. 

Thus  far,  we  have  talked  about  existence  results  of 
local  type.  This  means  that  we  have  initial  data  given  at  a 
certain  time,  and  we  prove  that  solutions  exist  on  some  time 
interval  thereafter,  which  may  be  short.  For  large  times, 
smooth  solutions  may  not  exist  as  we  shall  see  in  the  next 
section.  However,  there  are  results  on  existence  of  smooth 
solutions  for  all  times  if  one  considers  smooth  and  small 
data  perturbing  a  stable  state  of  rest.  Dafermos  ana  Nohel 
[5],  [6]  have  established  such  a  result  for  the  problem 


equation  will  prop, 
nonlinearities  are 
lop  from  initially 
iediate  nature  o£  mi 
eads  to  interestin' 


In  recent  joint  work  with  W.  Hrusa  [161,  we  study  (2.1) 
for  x  €  R  with  the  following  initial  data 

u  *  0,  t  <  0 

u(x,t  =0)  =  ^  o  sgn  x  ,  (2.2) 

ut(x,t  =  0)  =  6  sgn  x  . 

The  problem  (2.1),  (2.2)  can  easily  be  solved  by  Laplace 
transforms.  We  are  interested  in  the  regularity  of  the 
solution  and  the  location  of  singularities.  It  is  easy  to 
see  that  these  are  determined  by  the  asymptotic  behavior  at 
infinity  of  the  Laplace  transform  of  u,  which  is  related 
to  the  transform  of  m. 

If  m  is  a  smooth  function,  then,  for  A 
Re  X  >  0,  the  transform  has  the  asymptotic  behavior 

m(  A )  «  i  m( 0)  +  0(±~)  .  (2.3) 


In  this  case,  it  is  well  known  [1]  that  the  singularity 
propagates  with  the  wave  speed  /b  +  m ( 0 )  and  decays 

A 

exponentially  with  a  factor  exp(-m( 0)t/2(b  +  m(0))).  In 
addition,  we  show  in  116]  that,  except  in  a  very  special 
case,  there  is  a  discontinuity  in  uxx  across  x  *  0.  Away 

from  the  lines  x  =  ±t/b  +  m( 0) ,  x  *  0,  the  solution  is 
smooth. 

The  presence  of  a  singularity  in  m  alters  the 

A 

asymptotic  behavior  of  m(\).  Of  particular  physical 
interest  are  kernels  m  which  are  smooth  away  from  0,  but 
become  infinite  at  0.  Such  a  possibility  has  been 
suggested  by  certain  molecular  theories  [7],  [34],  [39]  and 
by  experiments  [22]. 

As  is  expected  from  the  formula  for  the  wave  speed 
above,  there  remains  a  finite  wave  speed  as  long  as  m  is 
integrable.  However,  the  exponential  decay  factor  becomes 
zero  if  m(0)  is  infinite,  and  we  can  no  longer  expect 


propagation  of  a  discontinuity,  in  fact/  Hannsgen  and 
Wheeler  110]  show  that,  for  any  singular  kernel,  there  must 
be  some  kind  of  "smoothing";  that  is,  the  solution  at  any 
positive  time  will  lie  in  a  more  restricted  function  space 
than  the  space  admitted  for  the  initial  data.  How  much 
smoothing  there  is,  depends,  not  surprisingly,  on  the  nature 
of  the  singularity  in  m.  In  [16] ,  [28] ,  the  following 
examples  are  discussed: 


-n°t 


1.  m(t )  =  £  6  9  a>w# 

n=i  , 

This  kernel  behaves  like  t  /<x  as  t  ♦  0,  and  is 
integrable  if  a  >  1.  The  solution  is  C-  across  the  wave 


and  analytic  behind  the  wave,  i.e.  for  0  <  x  <  t  /b  +  m(0) 
(the  latter  is  in  fact  true  for  any  completely  monotone 
kernel).  If  a  <  1,  m  is  not  integrable.  The  wave  speed 
is  infinite,  and  the  solution  is  analytic  everywhere  in 
t  >  0,  x  *  0. 


2.  m(t)  =  l  exp( -ent )  . 
n=l 

This  kernel  behaves  like 


in  t 


as  t  0. 


number  of  derivatives  of  the  solution  which  exist  across  the 
wave  is  increasing  proportional  with  time. 


3.  m(t)  *  /  exp(-ee  t)dv  . 

0 

This  behaves  like  in  I  in  t 


This  behaves  like  in | in  t|  as  t  ♦  0.  If  regularity 
is  measured  by  the  number  of  existing  derivatives,  the  gain 
of  regularity  in  the  solution  is  infinitesimal.  However, 
the  initial  discontinuity  disappears  and  u  is  continuous 
across  the  wave. 

We  believe  that  continuity  across  the  wave  should  hold 
under  very  mild  assumptions  on  m  (provided  of  course 
that  m  is  singular),  but  we  have  not  been  able  to  prove  a 
general  result.  It  would  be  enough  to  show  that  the  total 
variation  remains  finite. 

As  this  discussion  shows,  singular  kernels  lead  to 
completely  new  qualitative  features  of  wave  propagation. 

The  coexistence  of  finite  speed  and  smoothing  may  be  of 


VS 


interest  in  modeling  phenomena  in  other  fields  of  physics, 
such  as  heat  flow.  It  would  be  extremely  interesting  to 
know  how  singular  kernels  affect  the  solutions  of  nonlinear 
problems,  but  no  results  exist  so  far.  In  the  following,  I 
shall  review  some  results  for  nonlinear  problems  with  smooth 
kernels. 

Coleman  and  Gurtin  [3],  partly  in  joint  work  with 
Herrera  [2]  have  studied  the  evolution  of  acceleration 
waves.  They  find  that,  if  an  acceleration  wave  propagates 
into  a  medium  at  rest,  then  its  amplitude  satisfies  a 
Riccati  equation  of  the  form 

gf  =  _Ya  +  ^  a2  .  (2.4) 

Here,  in  our  notation  from  above,  y  *  - m-  — i - ,  and  A 

2 ( b  +  m( 0) ) 

is  a  second  constant  related  to  nonlinear  properties  of  the 
material.  From  (2.4)  it  is  clear  that  a  will  decay 
exponentially  if  its  initial  value  is  small,  but  will  blow 
up  in  finite  time  if  its  initial  value  is  large  and  has  the 
right  sign. 

This  result  shows  that  weak  singularities  can  grow  in 
amplitude  and  lead  to  stronger  singularities.  This,  and 
known  results  on  hyperbolic  conservation  laws,  have 
motivated  the  search  tor  results  establishing  the 
development  of  singularities  from  smooth  data.  For  certain 
specific  model  equations,  such  results  were  obtained  by 
Slemrod  135],  Gripenberg  (8],  Hattori  19],  Malek-Madani  and 
Nohel  [24]  and  Malek-Madani,  NOhel  and  John  [25].  We  do  not 
present  these  proofs  here,  but  the  following,  somewhat 
sketchy  argument  will  illustrate  some  of  the  main  ideas. 

Let  us  consider, the  model  problem  studied  by  Malek-Madani 
and  Nohel  [24] 

t 

ut  +  4(u)x  -  /  m(t  -  s)t(u(s))xds  =  0  .  (2.5) 

Here  4  and  4  are  strictly  monotone  functions: 


-8- 


»  •  «  «,  » 


.  *S.  -V 


♦',♦'>  0  It  we  set 


2  =  -  /  m(t  -  s)iMu(s)  )ds  , 
0 

then 


z¥  =  -m(0)*(u)  -  /  m'(t  -  s)\|i(u(s)  )ds 
c  0 


(2.7) 


If  t  is  small  and  u  is  bounded,  then  the  integral  in 
(2.7)  is  small.  We  can  then  regard  (2.5)  as  a  perturbation 
of  the  problem 


ut  +  ♦(u)x  +  zx  =  0 
zt  =  -m(O)iKu)  , 


(2.8) 


and  we  expect  that  the  solution  of  (2.5)  will  develop  a 
discontinuity  in  u  if  the  solution  to  (2.8)  develops  such 
a  discontinuity  sufficiently  fast,  (and  in  fact,  it  is  not 
difficult  to  include  the  integral  as  a  perturbation  in  the 
following  argument).  In  (2.8),  we  introduce  Riemann 
invariants  r«t(u)+z,  s  *  z,  and  we  define  the 
characteristic  derivative  +  ♦ '  (u)  |^  .  We  then 

obtain 


g£  -  -m(0)*(u)  , 


||  =  -m(O)iMu)  . 


(2.9) 


We  differentiate  with  respect  to  x,  and  use  the  fact  that 
$(u)  =  r  -  s,  to  obtain 


Sr  r  .  -  t*ig>  (r  -  sir  -  »(Q)  »'<“>.  (r  -  s) 
dt  x  $  * (u)  x  xx  ♦ '  ( u )  x  x 


TE  sx  =  _m(0)  f*TiTT  (rx  "  sx} 


(2.10) 


If  we  now  let  p  =  sup  |r  |,  and  a  =  sup  J s  | ,  then, 

x  x  A 

under  appropriate  hypotheses  on  $,  y  and  the  initial  data 
(2.10)  leads  to  differential  inequalities  of  the  form 


v-VVVWV 


(gYJ+p  >  ap  ~  ®p0  “  Y(p  +  ®) 
(|^)+0  <  6(P  +  o) 


(2.11) 


It  is  not  hard  to  show  that  we  can  take  data  such  that  p 
tends  to  *  in  a  short  time. 

The  analogy  with  hyperbolic  conservation  laws  and  the 
numerical  evidence  [26]  suggest  that  this  blow-up  of  smooth 
solutions  means  the  development  of  shocks  (i.e.  weak 
solution  with  discontinuous  u  for  (2.5)  or  solutions  with 
discontinuous  ux  and  ut  for  second  order  equations).  N 
such  result  has  been  proved.  Existence  theorems  for 
solutions  with  initial  data  of  small  total  variation  were 
established  by  Dafermos  and  Hsiao,  but  their  work  is  still 
unpublished . 

3.  STEADY  FLOWS  OF  VISCOELASTIC  FLUIDS 


Existence  results  for  steady  flows  of  viscoelastic 
liquids  are  lacking,  even  for  slow  flows  perturbing  rest. 
Formally,  such  slow  flows  have  been  analyzed  by  perturbation 
methods.  Let  the  stress  be  given  by  a  smooth  functional  of 
the  history  of  the  relative  Cauchy  strain 


t 

T  =  F 
s=-° 


lCt(s))  . 


(3.1) 


We  assume  Ct(s)  is  a  smooth  function  of  s,  expandable  in 
a  Taylor  series  about  s  =  t: 


£*(8)  -  1  +  l 
c  i  =  l 


(s  -  t)' 


^(t)  +  0(  |s  -  t 


)  .  (3.2) 


The  Rivlin-Ericksen  tensors  satisfy  the  following 

recursion  relation  in  a  steady  flow 

A .  ,  »  ( u* V ) A.  +  A.Vu  +  (Vu)TA.  , 

-l+l  -  _i  -i  -  -  -l  ( 3.3) 

=  Vu  +  (Vu)T  . 

We  can  therefore  regard  (3.2)  as  an  expansion  of 
£fc(s)  in  increasing  powers  of  a  small  and  smooth  velocity 
field  u.  Suppose,  e.g.,  that  we  want  to  solve  the 


-10- 


equations  of  motion 


p(u*V)u  »  div  J  -  7p  +  ef  , 

div  u  =  0  ,  (3.4) 


in  a  bounded  domain  fl,  where  e  is  a  small  parameter. 
Under  appropriate  conditions  on  the  functional  F  (see 
Coleman  and  Noll  [4]),  we  can  formally  expand  u  in  powers 
of  e  and  equate  coefficients  of  the  same  powers  of  e  in 
the  equations.  This  leads  to  an  infinite  sequence  of 
inhomogeneous  Stokes  problems  from  which  an  approximate 
"solution"  can  be  determined  to  any  order  in  e. 

It  has  not  been  proved  that  such  a  formal  expansion 
converges  or  that  it  yields  an  asymptotic  approximation  to  a 
solution  of  the  full  equation.  Because  of  the  presence  of 
the  term  (u*V)  ^  in  (3.3),  the  expansion  of  Qt(s) 
involves  derivatives  of  u  of  arbitrary  order.  It  is 
therefore  clear  that  we  can  expect  convergence  of  the 
expansion  only  if  we  require  a  high  degree  of  smoothness  for 
the  data  of  the  problem.  It  is  possible  that,  under 
appropriate  assumptions,  convergence  could  be  shown  in 
certain  spaces  of  analytic  functions.  An  effort  in  this 
direction  was  made  by  Niggemann  [27].  He  studies  one- 
dimensional  model  problems  of  the  form 

-  -  (ji)  <3n> 

u"  *  f  +  l  l  a,  .  uu  ...  u  (3.5) 

"-1  h . v1  1"  N 

His  analysis  requires  certain  assumptions  on  the 

coefficients  a  .  ,  for  which  no  physical  motivation 
D  1*  •  On 

has  been  given. 

An  alternative  approach  to  the  problem  is  the  use  of 
iterative  schemes,  by  which  it  may  be  possible  to  prove 
existence  of  solutions  without  recourse  to  small  parameter 
expansions.  The  asymptotic  validity  of  such  expansions  can 
then  be  establihsed  a  posteriori,  by  estimating  derivatives 


of  the  solution  with  respect  to  e.  In  recent  work  by  the 
author  [33]#  existence  of  solutions  to  (3.4)  is  proved  in 
this  way,  provided  that  uQ  is  tangential  to  3Q#  and  j 
is  determined  by  a  differential  constitutive  law  of  Maxwell 
or  Jeffreys  type#  with  an  arbitrary  number  of  relaxation 
modes.  For  simplicity#  we  hall  here  consider  only  the 
special  case  of  the  upper-convected  Maxwell  model 


(U-V)|  -  (Vu)J  -  |(7u)T  +  X|  =  nX(7u  +  (Vu)T) 


(3.6) 


By  taking  the  divergence  of  this  equation#  we  obtain 
(u«V)div  |  -  (Vu)div  T  +  X  div  T  *  T  :  32u  +  nXAu  #  (3.7) 

2 


2  3 

where  T  :  3  stands  for  V  T ■ .  ... ■  .  We  substitute 

j,k  3  sVxk 


div  X  =  p(u*V)u  +  Vp  -  f  from  (3.4)  (we  set  e  =  1), 
introduce  the  notation  q  =  (u»7)p  +  Xp,  and  solve  the 
equations  by  the  following  iteration  scheme. 


y0  -  Q,  P°  -  q°  =  0#  X°  -  0 


:0 

»n 


(3.8) 


T"  s  32un+1  +  nXAun+1  -  p(un«V) (un»V)un+1  -  Vqn+1 


=  -  [Vun  +  (Vun)TJVpn  -  ( u"  •  V )  f  +  (Vun)f  -  Xf 


(3.9) 


-p(?un 

)(un«V)un  +  Xp (un»  V ) un 

div  un+1  =  0, 

Sn+1(sa  ’  H0'  ///  o"+1  ■  « 

( un+3-  •  V  )pn+^  + 

xpn+l  .  qn+l 

(3.10) 

(un+1-V)Xn+1  - 

+  XJn+1  = 

(Vun+1)Xn+1  -  |n+1(Vun+1)T 

nX(Vun+1  +  (Vun+1)T) 

(3.11) 

Equation  (3.9)  is  basically  a  perturbation  of  the  Stokes 
problem#  while  (3.10)  and  (3.11)  are  hyperbolic  problems 
which  can  be  solved  by  the  method  of  characteristics.  The 
proof  for  convergence  of  (3. 8) -(3.11)  is  similar  to 
existence  proofs  for  initial  value  problems  for  quasilinear 


hyperbolic  systems.  First  bounds  for  all  iterates  in  a 
sufficiently  strong  norm  are  obtained ,  and  then  these  bounds 
are  used  to  obtain  a  contraction  estimate  in  a  weaker 
norm.  Iterations  similar  to  (3.8)-(3.11)  have  been  used  in 
numerical  calculations  (see  e.g.  [40]). 

The  following  questions  remain  open: 

1.  Can  similar  results  be  obtained  for  integral  type 
constitutive  models? 

2.  What  happens  if  uQ  is  not  tangential  to  ao?  In 
this  case  we  expect  from  the  analysis  of  characteristics 
[18]  that  extra  boundary  conditions  need  to  be  imposed  at 
inflow  boundaries,  fthat  are  the  right  conditions  to  impose 
and  what  compatibility  conditions  are  required  to  avoid 
singularities  of  the  solution  where  an  inflow  boundary  joins 
with  other  parts  of  the  boundary? 

3.  Formal  expansions  for  small  parameters  require  no 
extra  boundary  conditions  at  inflow.  How  are  solutions 
approximated  by  these  expansions  distinguished  from  others? 

4.  Can  one  say  anything  about  steady  flows  when  f 
and  Uq  are  not  small?  In  this  case  it  is  possible  that  a 
change  of  type  occurs  in  the  equations  governing  steady  flow 
[18].  However,  attempts  to  compute  steady  flows  have 
encountered  difficulties  even  in  situations,  such  as 
creeping  flows  of  an  upper  convected  Maxwell  model,  where 
there  should  be  no  change  of  type.  Although  the  existence 
of  steady  flows  in  the  Newtonian  case  is  well-known  [38], 
even  at  high  Reynolds  number,  we  have  no  reason  to  believe 
that  steady  flows  of  viscoelastic  fluids  should  always 
exist.  It  is  possible  that  the  numerical  difficulties 
result  from  non-existence. 
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